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Aiming to understand the main spin relaxation mechanism in graphene, we investigate the spin 
relaxation with random Rashba field induced by both adatoms and substrate, by means of the 
kinetic spin Bloch equation approach. The charged adatoms on one hand enhance the Rashba spin- 
orbit coupling locally and on the other hand serve as Coulomb potential scatterers. Both effects 
contribute to spin relaxation limited by the D'yakonov-Perel' mechanism. In addition, the random 
Rashba field also causes spin relaxation by spin-flip scattering, manifesting itself as an Elliott- Yafet- 
like mechanism. Both mechanisms are sensitive to the correlation length of the random Rashba field, 
which may be affected by the environmental parameters such as electron density and temperature. 
By fitting and comparing the experiments from the Groningen group [Jozsa et al, Phys. Rev. 
B 80, 241403(R) (2009)] and Riverside group [Pi et al, Phys. Rev. Lett. 104, 187201 (2010); 
Han and Kawakami, ibid. 107, 047207 (2011)] which show either D'yakonov-Perel'- (with the spin 
relaxation rate being inversely proportional to the momentum scattering rate) or Elliott- Yafet-like 
(with the spin relaxation rate being proportional to the momentum scattering rate) properties, we 
suggest that the D'yakonov-Perel' mechanism dominates the spin relaxation in graphene. The latest 
experimental finding of a nonmonotonic dependence of spin relaxation time on diffusion coefficient 
by Jo et al. [Phys. Rev. B 84, 075453 (2011)] is also well reproduced by our model. 

PACS numbers: 72.25.Rb, 71.70.Ej, 67.30.hj, 05. 40. -a 



I. INTRODUCTION 



Graphene, a two-dimensional allotrope of carbon 
with a honeycomb lattice, has attracted much at- 
tention due to its two dimensionality, Dirac-like en- 
ergy spectrum and potential for the all-carbon based 
electronics and spintronics in recent years**- With 
the breaking of inversion symmetry, possibly caused 
by ripples, 9 perpendicular electric fields^r— adsorbed 
adatoms,— ~— the substrate,— ~— etc., the Rashba spin- 
orbit couplin g 10 ' 11 ' 19 arises and results in spin relaxation 
in the presence of scattering in graphene. A number of 
experiments on spin relaxation in graphene on SiO*2 sub- 
strate are available, revealing a spin relaxation time of 
the order of 10-100 ps^SSr— However, somewhat contra- 
dictory results are exhibited in these experiments. While 
a decrease of spin relaxation rate with increasing mo- 
mentum scattering rate has been observed by Riverside 
group by surface chemical doping at 18 K^i a linear 
scaling between the momentum and spin scattering has 
been observed by both Groningen group 20 ' 22 : 23 at room 
temperature and very recently Riverside group 2 ^ at low 
temperature (< 100 K) via varying the electron density 
in graphene. The D'yakonov-Perel— (DP) mechanism 
was justified to be important by the former phenomenon, 
however, the Elliott- Yafet^ 9 - (EY) mechanism was sug- 
gested to account for the latter. In addition, at a tem- 
perature as low as 4.2 K, a nonmonotonic dependence 
of spin relaxation time on diffusion coefficient with the 
increase of electron density was also reported by Jo et 
al. very recently^ They claimed that the spin relax- 



ation is due to the EY mechanism. To fully understand 
the spin relaxation in graphene, theoretical studies are in 
progress! 15 i 18 i 30 ~— 

Theoretically, the EY mechanism is revealed to be 
both invalid to account for the linear scaling between 
the momentum and spin scattering with the increase 
of electron density^ 3 - and also less important than the 
DP one. 30 According to Refs. HO and IH, the spin re- 
laxation times caused by the EY and DP mechanisms 
are rf Y w (hvfkf /X) 2 t p and r s DP « h 2 /(X 2 T p ) respec- 
tively. Here Vf and kf are the Fermi velocity and mo- 
mentum, A is the Rashba spin-orbit coupling strength 
and t p is the momentum relaxation time. Consequently, 
rf Y /t p oc fc 2 ) 30 ' 33 disagreeing with the electron density- 
independent linear scaling between t s and t v in the 
experiments! 20 ' 22 ' 23 ' 26 The latest experiment by Jo et al. 
reveals the relation rf Y /r p oc k 2 with the increase of 
carrier density and hence suggests that the EY mecha- 
nism dominates spin relaxation in graphenei^l Nonethe- 
less, from the above theoretical results it is found that 
rf Y /t,P p s» (kfl) 2 with I being the mean free path. Typi- 
cally (kfl) 2 ^ 1 in graphene, meaning that the DP mech- 
anism dominates over the EY one^°- Due to the above 
reasons, the EY mechanism can not dominate spin relax- 
ation in graphene. 

The DP mechanism then becomes the reasonable can- 
didate for the dominant spin relaxation mechanism. Ini- 
tially, the spin relaxation time limited by the DP mech- 
anism was calculated to be much longer than the ex- 
perimental data ~10-100 ps due to the weak Rashba 
field. For example, ripples with curvature radii ~ 50 nm 
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induce the local Rashba spin-orbit coupling with A ~ 
8.5 //eV, 9 and a perpendicular electric field of magnitude 
E z (E z ~ 0.1 V/nm) contributes a Rashba spin-orbit 
coupling with A ~ (E z , where C is 0.258 /j,eV/(V/nm) 
from rough estimation^ 17.9 or 66.6 ^eV/(V/nm) from 
the tight-binding model 9,11 and 5 /j,eV/(V/nm) from 
the first-principles calculation ! 12 ! 13 Based on this weak 
Rashba spin-orbit coupling, Zhou and Wu calculated spin 
relaxation in graphene on Si02 substrate with mobil- 
ity ~ 10 4 cm 2 /(V-s) taken from the charge transport 
measurement^ and obtained a quite long spin relaxation 
time of the order of fj,sM- However, for the nonlocal mea- 
surements of the spin relaxation, the mobility is at least 
one order of magnitude smaller,— ~ 23 > 26 i 27 most likely 
caused by the extrinsic factors induced by the ferromag- 
netic electrodes^ e.g., the adatoms. The adatoms, as 
well as the influence of the substrate, may substantially 
enhance the Rashba spin-orbit coupling by distorting the 
graphene lattice and inducing sp 3 hybridization, leading 
A to be ~ meVM r i 1 ' With this enhanced Rashba spin- 
orbit coupling, the spin relaxation time in graphene is 
estimated^ and calculated^ to be comparable to the 
experimental data. Nevertheless, whether and how the 
DP mechanism accounts for the experimentally observed 
linear scaling between the momentum and spin scattering 
is questionable. 

Apart from the above two mechanisms, another EY- 
like mechanism may contribute to the spin relaxation in 
graphene when the Rashba field is random in the real 
space. As proposed by Sherman in semiconductors ^L22. 
the randomness of spin-orbit coupling contributes to or 
even dominates spin relaxation by spin-flip scattering un- 
der certain conditions i 32 i 36 ~ ^S For graphene, the Rashba 
field induced by a fluctuating electric field from ionized 
impurities in the substrate or ripples is indeed random in 
the real space. The former case, with the average Rashba 
field being nonzero, has been investigated by Ertler et al. 
via Monte Carlo simulation,— while the latter one, with 
the average Rashba field being zero, has been studied by 
Dugaev et al£2- via the kinetic equations. However, for 
both cases the calculated spin relaxation time is much 
longer than the experimental data. 

In this work, we investigate spin relaxation in graphene 
with random Rashba field (RRF) induced by adatoms 
and substrate by means of the kinetic spin Bloch equa- 
tion (KSBE) approach^ A random Rashba model is set 
up, where the charged adatoms on one hand enhance 
the Rashba spin-orbit coupling locally and on the other 
hand serve as Coulomb potential scatterers. Based on 
this model, an analytical study on spin relaxation with 
RRF is performed. It is found that while the average 
Rashba field leads to spin relaxation limited by the DP 
mechanism, which is absent in the work by Dugaev et 
al.j^ the randomness causes spin relaxation via spin-flip 
scattering. With the increase of adatom density, the spin 
relaxation caused by the spin-flip scattering due to the 
RRF always shows an EY-like behaviour (the spin relax- 
ation rate is proportional to the momentum scattering 



rate) whereas the DP mechanism can exhibit either EY- 
or DP-like (the spin relaxation rate is inversely propor- 
tional to the momentum scattering rate) one. When all 
the other parameters are fixed, with the increase of elec- 
tron density the spin relaxation rates due to both mech- 
anisms increase; Nevertheless, the spin relaxation rate 
determined by the spin-flip scattering due to the RRF is 
insensitive to the temperature whereas that determined 
by the DP mechanism becomes insensitive to the tem- 
perature when the electron-impurity scattering is domi- 
nant. However, the correlation length of the RRF may 
vary with the electron density as well as temperature and 
both mechanisms are sensitive to the correlation length. 

We carry out numerical calculations and fit the ex- 
periments of Riversid o 24 ' 26 and Groningen 23 groups. By 
fitting the DP-like behaviour with the increase of adatom 
density observed by the Riverside group^ we find that 
only when the DP mechanism is dominant can the exper- 
imental data be understood and the effect of the spin- 
flip scattering due to the RRF is negligible. However, 
the experimental EY-like behaviour with the increase of 
electron density first observed by the Groningen group 23 
can be fitted from our model with either the DP mech- 
anism or the spin-flip scattering due to the RRF being 
dominant by taking into account the decrease of the cor- 
relation length of the RRF with the increase of electron 
density. Nevertheless, the fact that the Riverside group 
has also observed the similar EY-like behaviour in their 
samples very recently 2 - suggests that the DP mechanism 
is dominant but exhibits EY-like properties with the in- 
crease of electron density. The temperature dependence 
of the spin relaxation from the Riverside group,— sug- 
gested to be the evidence of the EY mechanism, is also 
fitted by our random Rashba model with the DP mech- 
anism being dominant. The corresponding temperature 
dependence from the spin-flip scattering due to the RRF 
is demonstrated to be in fact temperature insensitive. 
The similar experimental phenomenon observed with the 
variation of the electron density by the two groups fur- 
ther suggests that the DP mechanism also dominates the 
spin relaxation in the experiment of Groningen group^ 
Moreover, the latest reported nonmonotonic dependence 
of spin relaxation time on diffusion coefficient^ 7 - is also 
fitted. 

This paper is organized as follows. In Sec. II, we 
present the model and introduce the KSBEs. In Sec. Ill, 
we investigate the spin relaxation analytically and dis- 
cuss the relative importance of the DP mechanism and 
mechanism of the spin-flip scattering due to the RRF. 
In Sec. IV, we perform numerical calculations and fit the 
experimental data. We discuss and summarize in Sec. V. 



II. MODEL AND KSBES 

The n-doped graphene monolayer under investigation 
lies on the Si02 substrate perpendicular to the z-axis. 
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The random Rashba spin-orbit coupling read o 10 i n 

H so = \(r)((iT x a y - T y a x ). (1) 

Here // = ±1 labels the valley located at K or if', r 
and cr are the Pauli matrices in the sublattice and spin 
spaces, respectively. The position-dependent coupling 
strength A(r), mainly contributed by the randomly dis- 
tributed adatoms and also possibly by the substrate, can 
be modeled as 

N a 

A(r)=A*+5> ie -l r - R "l 2 / 2 « 2 . (2) 

n=l 

Here the second term is contributed by the adatoms 
with a total number Nf. In this model it is assumed 
that an adatom located at R„ induces a local Rashba 
field peaking at R„ with a magnitude of 8 n and de- 
caying within a length scale ~ £ following the Gaussian 
form. S n is of the order of meV (Refs. Il3l - [l5l) and £ is 
larger than the graphene lattice constant 0.25 nm^ The 
first term Aq comes from the average contribution from 
the substrate, whose fluctuation is phenomenally incor- 
porated by affecting £. The mean value of A(r) reads 
A = (A(r)) = A + Ag with Ag = 2tto?&; 2 , where nf is 
the areal density of adatoms. The correlation function 
C(r) = ([A(r) - A o ][A(0) - A ]> - imfS s ?e-'*'< a , with 
the corresponding Fourier transformation 

C7 q = 4^ 2 <^ 4 e^ V . (3) 

This expression is similar to that given in Ref. |32| ex- 
cept that it depends on £ in a higher order here. For 
single-sided adatoms (the adatoms are distributed on 
the graphene surface), we choose 8 n = 8, by which 
8 = 8 and S 2 = S 2 . For double-sided adatoms (the 
adatoms are distributed both on the graphene surface 
and at the graphene/substrate interface), we set 8 n = Si 
or 82 {8182 < 0) with equal possibilities, by which 
6 = 0.5(5i + 6 2 ) and P = 0.5(<5 2 + <5f). It is noted 
that the random Rashba model proposed here is modified 
from the short-range random potential model depicting 
the electron-hole puddles in graphene 41 

Under the basis laid out in Refs. andHH, the electron 
Hamiltonian can be written as^ 1 - 

H = ^ [tkSss' + J^k ' Owjc^ks^ks' + H in t (4) 
/^kss' 

in the momentum space. Here c^ks (c M ks^) is the anni- 
hilation (creation) operator of electrons in the p valley 
with momentum k (relative to the valley center) and spin 
s (s = ±i). £k = hvfk with vj = 10 6 m/s. The effective 
magnetic field J~ik from the average Rashba field is 

Ok = Ao(-sin0k,cos0k,O), (5) 

where (9k is the polar angle of momentum k. The 
Hamiltonian H lnt consists of the spin-conserving scat- 
tering [electron-impurity^ 2 - (here the impurities in- 
clude both the ones existing in the substrate and the 



charged adatoms ) i taken into account by the minimal 
model proposed by Adam and Das Sarmai 2 -), electron- 
phonon,— ~— and electron-electron 3 - 1 - scattering] as well 
as the spin- flip scattering due to the RRF t 32 i 38 

-ffflip = ^ Ak-k'Vkk'ss'C^ks^W' (6) 
/^,k'7^k,ss' 

where 

A q = J A(r)e- lq r dr (7) 

and 

= ( X ) ■ (8) 

The KSBEs are42 

9 t p M k(t) = <9tp M k(i)|coh + 5 t^k(i)lscat + d t p^{t)\{ cat . 

(9) 

Here ptfk{t) represent the density matrices of electrons 
with relative momentum k in valley p at time t. The 
coherent terms d t p^{t)\ coa = -|[O k ■ er, p^t)] with 
the Hartree-Fock term from the Coulomb interaction 
being neglected due to the small spin polarizationi 31 i 40 
The concrete expressions of the spin-conserving scatter- 
ing terms 9tP^k(*)|scat can be found in Ref. |3l|. When 
the electron mean free path I is much longer than the 
correlation length of the fluctuating Rashba field £ (this 
is easily satisfied as £ ~ nm while I ~ 10-100 nm in 
graphene), i.e., the electron spins experience indeed the 
random spin-orbit coupling, the spin-flip scattering terms 
can be written as (Appendix [A" f 32 i 38 

k' 

X [p M k(t)-Vkk'P M k'(*) F k'k]- (10) 

By solving the KSBEs, one can obtain the spin relaxation 
properties from the time evolution of p^\s.{t). 

III. ANALYTICAL STUDY OF SPIN 
RELAXATION 

In this section we analytically study the spin relaxation 
in graphene with the RRF. To realize this, we only take 
into account the spin-flip scattering as well as the strong 
elastic electron-impurity scattering. When the valley in- 
dex is further omitted due to the degeneracy, the KSBEs 
are simplified to be 

i 27T 

d t pu{t) = --[n k - - — ^|[/ k - k 'l 2 4k' 

k' 

x 6(e k - £k')[Pk(i) - pk'(t)\ 

--jr^2 Ck-k'£(£k - £k ,) 
k' 

x [pk(t)~ Vkk'/0k'(i)Vk' k ] (11) 
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Here |£/k-k'| 2 is the effective efectron- impurity scatter- 



ing matrix efement and / kk / = |[1 + cos ( 



0] is 



^k-k<l 2 > 



the form factor 41 |J7 k -k'| 2 = nf\U£_ k , 
with the first (second) term corresponding to the scat- 
tering of electrons from adatoms (impurities in the sub- 
strate). |C^k-k'l 2 i s t Qe electron-impurity Coulomb po- 
tential scattering matrix element i 31 i 34 ' 42 n\ is the impu- 
rity density in the substrate. By defining the spin vector 
as S k (i) = Tr[p k (i)er], one obtains the equation of S k (t) 
directly from Eq. ([TT]) as 



ftS k (t)= ^Ak-SkW-^^lC/^l 2 /, 



kk' 



From Eq. (|T5| . one notices that the spin relaxation 
rate (take T z as an example) consists of two parts, 
Tflip = 2/tP s , determined by the spin- flip scattering due 

to the RRF, and Tdp = 4A 2 ,t 1 : /ti 2 , determined by the av- 
erage Rashba field due to the DP mechanism. These two 
mechanisms contribute to spin relaxation independently. 
It is noted that T^p obtained here is consistent with that 
given in Ref. [32| except the correlation functions are dif- 
ferent and Tdp is the one previously given in Ref. HU- In 
the following we discuss the spin relaxation due to the 
two mechanisms respectively and compare their relative 
importance. 



X S(£k-£k')[Sk(t)-S k ,(i)] 
~ T" ^2 c k-k>S(ek - £kO 

k' 

x [S k (t)-B kk , -S k -(i)]. 



Here 



A k = 






cos # k 



cos # k 
sin # k 
sin # k 



(12) 



(13) 



and 



B 



kk' 



COS( 

- sinf 



fk'J 
9 k ') 



— sm(W k 
cos(0 k -f 



.(14) 







By expanding S k (f) as S k (f) = £ ( S^.(i)e 4 ' 9k and re- 
taining the lowest three orders of S l k (t) (i.e., terms with 
I = 0, ±1), one obtains a group of differential equations 
of S A , 1 '°(i) (Appendix IB]). With initial conditions, these 
equations can be solved and the information on spin re- 
laxation is obtained from S®(t). We label the spin re- 
laxation rate along the x-, y-, or z-axis for states with 
momentum k as T x , T y or T z , respectively. One has 
(Appendix [5]) 



2I\ 



2E, 



where 



and 



_L_ _ k_ 



|C/ q | 2 sin 2 ' 



'fcs 



2irh 2 v f J 



with \U a \ 2 and C q depending only on |q| 



(15) 



(16) 



(17) 



2fcsin§. For 



a highly degenerate electron system in graphene, the spin 
relaxation is contributed by the spin-polarized electrons 
around the Fermi circle. Therefore, one can approxi- 
mately obtain the spin relaxation rate of the whole elec- 
tron system by replacing k with kf in Eq. (|15[) . 



A. Spin relaxation caused by the spin-flip 
scattering due to the RRF 

By utilizing Eq. Q , one can obtain the spin relaxation 
rate along the z-axis due to the spin-flip scattering as 

k f 2rr 
rmp = -rj— / d0C q 
™ v f Jo 



h 2 v f 

d,TT 2 nfP^ 



i e - 2 « 2fc // (2^ 2 ^) 

S^nfPk^ 
F ^kf) = £—f-G{Zk f )(18) 



h 2 vf h 2v f 

In the limits /c^^ <gc 1 and kf£ ^> 1, 



Tflip 



8w 2 n?8 2 



h 2 v f U 3 /(20F) 



(k f i « 1) 

{kfi » i) 



(19) 



In the above equations 1$ (x) is the modified Bessel func- 
tion, F(x) = xe~ 2x2 1 (2x 2 ) and G(x) = x 3 F(x). 

We now focus on the various factors affecting Tfljp. 
Tflip is proportional to the adatom density nf, as ex- 
pected. From the x dependence of F(x) and G(x) as 
shown in Fig. [TJ one can explore the dependence of T^p 
on kf and £, respectively. It is found that with the in- 
crease of kf, Tflip first increases almost linearly when 
kf£ < 0.63, then decreases mildly and eventually sat- 
urates [F(+oo) = -^j=\- As a result, with the increase of 
electron density n e , Tfn p has a nonmonotonic behaviour 
with a peak located at n c e = ^(^j^) 2 - When n e <C n% 

1/2 

Tflip is proportional to n e ', and when 



> n%, r flip be 



ip 

comes insensitive to n e . These features are in consistence 
with those presented in Ref. We give an estimation 
on n c e here. When the correlation length £ is set as 1 nm, 
rig w 1.3 x 10 13 cm -2 , which is a quite high value com- 
pared to the experimental data. Usually n e varies around 
10 12 cm~ 2 . Therefore, in order to observe the nonmono- 
tonic behaviour of Tfljp with increasing n e , £ is required 
to be around a relatively large value, e.g., 3.5 nm. When 
the electron density is fixed, Ffl ip increases monotoni- 
cally with increasing £, as indicated by the x dependence 
of G(x) in Fig. [1] The effect of temperature T on Tfljp 
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can be inferred from the k dependence of Tfu p . For the 
highly degenerate electron system, rg; p is insensitive to 
T. When the electron density is relatively low, with in- 
creasing T, Tflip is expected to increase mildly as elec- 
trons tend to occupy the states with larger momentum. 
With typical values S 2 = 25 meV 2 , n° = 3x 10 12 cm' 2 , 
n e = 10 12 cm -2 and £ = 0.5 nm, one has « 0.09 and 
Fgjp s» 670 ps. If £ is changed to be 2 times larger, i.e., 
1 nm, Tjjjp becomes about 16 times smaller, as Tflip is 
proportional to £ 4 when fc/£ <C 1 [Eq. (fl9l) ]. 




FIG. 1: (Color online) F(as) and G(x) in Eq. (T5J. The scale 
of G(x) is on the right-hand side of the frame. 



B. Spin relaxation caused by the DP mechanism 

In the analytical study, only the elastic spin-conserving 
scattering is considered. With the other spin-conserving 
scattering included, the spin relaxation rate due to the 
DP mechanism should be modified to be 



DP 



4A 2 r;(fc / )/fi 2 , 



(20) 



where T*{kf) is the effective momentum relaxation time 
limited by all the different kinds of scattering, includ- 
ing the electron-electron Coulomb scattering*^ When 
Ao is fixed, the electron density and temperature affect 
Tdp via Tp(kf). It has been shown previously 3 -^ that 
with the increase of n e or the decrease of T, the scat- 
tering strength decreases and Tdp increases. However, 
when the electron-impurity scattering is dominant, Tdp 
is insensitive to T. The dependence of Tdp on n° is 
not obvious. To facilitate the investigation, we write 

T p~ 1 ( k f) = T pJ( k f) + T p,a( k f), where T p,a( k f) is con- 
tributed by the electron-charged adatom scattering and 

T~i(kf) by all the other scattering. From Eq. ([To), one 

has 



T p,a( k f) = 



nfkf 



4:irh 2 Vf J 



with |q| = 2fc/sinf. With A = Aj + 2TrnfS£ 2 
c 2 n", Tdp can be written as 

4 (\h + c 2 n?) 2 



DP 



n2 r-l(k f ) + c in f 



(22) 



|[/ q a | 2 sin 2 , 



an. 



(21) 



which indicates a complex dependence on n". When 
c 2 = and A* ^ 0, c 2 > and < nf < c^Aj, - 
2c I lr p,U k f) or c 2 < and < < < -c^Aj, r DP 
decreases with increasing nf, exhibiting the DP-likc be- 
haviour. However, most interestingly, when c 2 > and 
nf > c^Xl - 2c^yl(k f ) or c 2 < and < > -^Xl 
Tdp increases with increasing nf, exhibiting the EY-like 
behaviour. Particularly, here in the limit with nf being 
large enough, Tdp <x n" approximately. With typical 
values in the presence of adatoms, 34 Ao = 0.2 meV and 
T*(kf) — 0.1 ps, r^p is calculated to be about 100 ps. 



C. Comparison of relaxations caused by the 
spin-flip scattering due to the RRF and the DP 
mechanism 

In this subsection we discuss the relative importance of 
the mechanism of spin-flip scattering due to the RRF and 
the DP mechanism in the regime with -C 1, which 
is typically realized in graphene. Under this condition, 
from Eqs. (19]) and (22|), one obtains r DP /r flip « lOtf/J 2 
when A° = T~l(k f ) = 0. For the case with single-sided 

adatoms, 5 = <5 2 , therefore Fop/Fflip 10 and Tflip can 
be neglected. However, for the case with double-sided 
adatoms, Tflip can be comparable to or even surpass T^p 
as S may be as small as zero. In reality T~l(kf) ^ and 
Tdp decreases with increasing T~j(kf). When the sub- 
strate also contributes to the Rashba field as A°, Tdp can 
be either enhanced (e.g., when \®5 > 0) or suppressed. 
As a consequence, for the configuration with single-sided 
adatoms, when the contribution from the substrate to 
the average Rashba field does not compensate that from 
the adatoms (e.g., when \®8 > 0) and the scattering 
other than the electron-adatom type is not extraordi- 
narily strong (i.e., t~\ is not unusually large), the spin 
relaxation caused by the spin-flip scattering due to the 
RRF can be neglected. In such a case, the spin relax- 
ation is limited by the DP mechanism with the adatoms 
contributing to the average Rashba field. This is just 
how the effect of adatoms was incorporated in our pre- 
vious investigation^ For other cases, whether the spin- 
flip scattering due to the RRF is important or not when 
compared to the DP mechanism is condition-dependent. 
Undoubtly, when the average Rashba field approaches 
zero, the spin-flip scattering due to the RRF tends to be 
important. In the work of Dugaev et al.^ the average 
Rashba field induced by ripples is zero and the spin re- 
laxation is solely determined by the spin-flip scattering 
due to the RRF. However, the spin relaxation time cal- 
culated in their model is of the order of 10 ns, two orders 
of magnitude larger than the experimental values. 
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IV. NUMERICAL RESULTS 

The KSBEs need to be solved numerically in order to 
take full account of all the different kinds of scattering. 
The initial conditions are set as 



n? (10 12 cm" 2 ) 



/Wo) 



fO , fO 
Jkf ^ Jki 



fO 
■I] 



f° 
J] 



kl - 

— n 



Tr[p Mk (0)n • <t] = n e P(0), 



//k 

53 Tr [^k(0)] = n e . 

//k 



At time t = 0, the electrons are polarized along n with 
the density and spin polarization being n e and P(fi), re- 
spectively, fkf/i * s the F errm distribution function of 
electrons with spins parallel/antiparallel to n, where the 
chemical potential is determined by Eqs. (I2"4"l) - (j2"5)) . By 
solving the KSBEs, one can obtain the time evolution of 
spin polarization along n as P{t) = Tr[p M k(i)n • 

<t] and hence the spin relaxation time t s . In the calcu- 
lation, we set .P(O) to be as small as 0.05 and n in the 
graphene plane, such as, along the y-axis, in order to 
compare with experiments. 



A. Adatom density dependence of spin relaxation 

In this section we study the adatom density depen- 
dence of spin relaxation based on the single-sided adatom 
model. In Fig.[5J the in-plane spin relaxation rate against 
the adatom density (at the top of the frame) or the in- 
verse of charge diffusion coefficient (at the bottom of the 
frame) is shown. The temperature is T = 300 K, the elec- 
tron density is n e = 10 12 cm~ 2 , the density of impurities 
in the substrate is n\ = 0.2 x 10 12 cm -2 and the param- 
eters for the single-sided adatom model are 8 — 5 meV 
and £ = 0.5 nm. In the figure, the spin relaxation rates 
with different values of Aq are plotted by the curves. The 
nearby data points of each curve are calculated with the 
spin-flip scattering being removed. The small discrep- 
ancy between each curve and corresponding data points 
indicates that the DP mechanism dominates the spin re- 
laxation. It is noted that in a large parameter regime 
of the background Rashba field A°, the curves show ob- 
vious EY-like behaviour, i.e., the spin relaxation rate is 
proportional to the momentum relaxation rate. When A° 
is large enough [larger than 2c2C^ 1 T~ I 1 (fc/) « 0.05 meV 
from the discussion in Sec lIIIB] , the spin relaxation rate 
decreases with increasing adatom density at low doping 
density of adatoms, showing the DP-like behaviour. 

We further apply the single-sided adatom model to 
reinvestigate the experiment of Pi et al. from Riverside^ 
which shows an obvious DP-like behaviour. At T = 
18 K, with the increasing density of adatoms (Au atoms) 
from surface deposition (although Au atoms also denote 
electrons to graphene, the electron density is fixed at 
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FIG. 2: (Color online) The dependence of in-plane spin relax- 
ation rate on the adatom density (at the top of the frame) or 
the inverse of the charge diffusion coefficient (at the bottom 
of the frame). The temperature is T — 300 K, the eiectron 
density is n e = 10 12 cm -2 , the density of impurities in the 
substrate is nf = 0.2 x 10 12 cm" 2 and the parameters for the 
singie-sided adatom model are S — 5 meV and £ = 0.5 nm. 
The curves are calculated with different Aq. For each curve, 
its nearby data points are calculated with the same parame- 
ters but without the spin-flip scattering. 



2.9xl0 12 cm~ 2 by adjusting the gate voltage 2 ^), the dif- 
fusion coefficient decreases while the spin relaxation time 
increases, as indicated by the crosses with error bars in 
Fig. [3l By fitting the experimental data without adatoms 
(before Au deposition), we obtain A° = 0.153 meV and 
n\ = 2.1 x 10 12 cm~ 2 . 34 During Au deposition, a group of 
parameters, S = 2.03 meV and £ = 0.5 nm, can reproduce 
the experimental data except when the adatom density 
is larger than 2xl0 12 cm -2 (the solid curve). By as- 
suming that £ decreases with the increase of nf when 
nf is large enough (the inset of Fig. [3]), the recalcula- 
tion can cover the experimental data in the region with 
large nf (the dashed curve). Similar to Fig. [21 the dots 
nearby the solid curve are calculated with the spin-flip 
scattering being removed (with £ being fixed as 0.5 nm). 
The small discrepancy between the solid curve and dots 
indeed indicates that the spin-flip scattering due to the 
RRF is not important. 



B. Electron density dependence of spin relaxation 

The electron density dependence of spin relaxation is 
studied by fitting the experiment of Jozsa et al. from 
Groningen, 23 which shows an EY-like behaviour. At 
room temperature, with the increase of electron density 
from 0.16 to 2.81 xlO 12 cm~ 2 (adjusted by the gate volt- 
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dence of diffusion coefficient, the impurity density in the 
substrate nf has to increase with increasing n e possibly 
due to the increased ionization (otherwise if nf is fixed, 
D will increase with increasing n e much more quickly), 
as shown by the dots in the inset. Meanwhile, with the 
increase of n e , £ should decrease as shown by the trian- 
gles in the inset (the scale is on the right-hand side of the 
frame) to account for the increase of r s . Otherwise if £ is 
fixed, t s will decrease with increasing D mainly due to the 
increase of Tdp, as indicated by the dashed curve in the 
figure. In Fig. O we also present a feasible fitting by the 
double-sided adatom model. The squares stand for the 
experimental data and the open circles are from calcula- 
tion. In our computation, Aq = 0, nf — 0.5 x I0 12 cm~ 2 
and Si = —S2 = 5 meV. The inset shows the dependences 
of £ (open triangles with the scale on the right-hand side 
of the frame) and nf (solid circles) on D when n e is in- 
creased. In this fitting, only the spin-flip scattering due 
to the RRF plays a role in spin relaxation. 



FIG. 3: (Color online) The dependence of in-plane spin relax- 
ation rate on the adatom density (at the top of the frame) or 
the inverse of the charge diffusion coefficient (at the bottom 
of the frame). The temperature is T — 18 K, the electron 
density is n e = 2.9 x 10 12 cm -2 and the density of impurities 
in the substrate is nf = 2.1 x 10 12 cm" 2 . The crosses with er- 
ror bars are the experimental data from Pi et al.^ The solid 
curve stands for the fitting data via the single-sided adatom 
model with Aj = 0.153 meV, S = 2.03 meV and £ = 0.5 nm. 
The nearby dots of the solid curve are calculated with the 
spin-flip scattering terms removed. The dashed curve is cal- 
culated with the same parameters as that of the solid curve 
except that £ decreases with nf when nf > 2 x 10 12 cm -2 , as 
shown in the inset. 



age) , both the charge diffusion coefficient and the spin re- 
laxation time increase, with the latter being proportional 
to the former (the squares in Fig. [4j . It should be noted 
that this EY-like behaviour can not be explained by the 
nearly linear curves shown in Fig. [2J as the linearity there 
is due to the increase of adatom density when the elec- 
tron density is fixed. In fact, according to Sec. IIII1 with 
the increase of n e as well as the accompanying increase 
of D [and hence the increase of T*(fe/)], both Tflip and 
Tdp should increase when the parameters for the adatom 
model are fixed [refer to Eqs. ([T9| and ((22)) ]. However, as 
will be shown in the following, with the assumption that 
£ decreases with increasing n e , both the single-sided and 
double-sided adatom models (hence both the DP mech- 
anism and the mechanism of the spin- flip scattering due 
to the RRF) are able to fit the experimental data. 

In Fig. 21 we present the fitting to the experimen- 
tal data of Jozsa et al. via the single-sided adatom 
model where the DP mechanism is dominant. The cal- 
culated spin relaxation time is plotted by the open cir- 
cles in Fig. |U The fitting parameters are chosen as 
Aq = 0.127 meV, nf = 0.5 x 10 12 cm" 2 and 5 = 4 meV 
here. In order to reproduce the electron density depen- 
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FIG. 4: (Color online) Fit to the experimental data (the de- 
pendence of spin relaxation time t s on the charge diffusion 
coefficient D with the increase of n e ) of Jozsa et aZ.— via 
the single-sided adatom model. The squares stand for the 
experimental data and the open circles are from our calcula- 
tion. The inset shows the dependences of £ (with the scale on 
the right-hand side of the frame) and nf on D when n e is in- 
creased. For comparison, the dashed curve is calculated with £ 
being fixed at 1.05 nm. nf = 0.5 xlO 12 cm -2 , \\ = 0.127 meV 
and 8 — 4 meV. 



C. Temperature dependence of spin relaxation 

We investigate the temperature dependence of spin re- 
laxation in graphene in this section. Although the spin 
relaxation time determined by the DP mechanism in- 
creases with growing T as pointed out in Sec. IIII1 this de- 
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FIG. 5: (Color online) Fit to the experimental data (the de- 
pendence of spin relaxation time r s on the charge diffusion 
coefficient D with the increase of n e ) of Jozsa et aZ.— via 
the double-sided adatom model. The squares stand for the 
experimental data and the open circles are from our calcula- 
tion. The inset shows the dependences of £ (with the scale on 
the right-hand side of the frame) and n- on D when n e is in- 
creased, nl = 0.5 xlO 12 cm -2 , Aj = and <5i = ~S 2 = 5 meV. 



pendence becomes very week when the electron-impurity 
scattering is dominant (which is satisfied in graphene on 
SiC>2 substrate), as revealed in Fig. 1 of Ref. |3l| (note the 
mobility there is even one order of magnitude larger than 
the ones in this investigation). The spin relaxation time 
determined by the spin-flip scattering due to the RRF 
is also insensitive to T, as shown in Fig. [6j Therefore, 
when other parameters are fixed, the spin relaxation in 
graphene is expected to depend on temperature weakly. 

It is quite interesting that a decrease of r s with T is 
observed by the Riverside group very recently^ More- 
over, when T is fixed, with the increase of n e (adjusted 
by the gate voltage), both D and t s increase, similar to 
the observations by Jozsa et al.£2- The decrease of r s 
with growing T may be due to the increase of the corre- 
lation length £ of the RRF with the increase of T, with 
either the DP mechanism or the spin-flip scattering due 
to the RRF being dominant. The linear scaling between 
t s and D with the variation of electron density also can 
not determine which mechanism is the dominant one, as 
demonstrated in the previous section. 

As a feasible way, we fit the temperature dependence of 
spin relaxation based on the single-sided adatom model 
by assuming £ to increase with T. One possible fit- 
ting is shown in Fig. [7J Experimentally, when the 
gate voltage Vcnp = 20 V (60 V), the electron density 
n e = 1.47 x 10 12 cm- 2 (4.42 x 10 12 cm" 2 ). 47 In Fig.^a) 
and (b), the squares (open circles) are the experimental 
data of spin relaxation time and diffusion coefficient with 



FIG. 6: (Color online) Temperature dependence of spin relax- 
ation time due to the spin-flip scattering. The parameters are 
n e = 1.5 x 10 12 cm' 2 , n° = 0.5 x 10 12 cm" 2 and S 2 = 4 meV 2 . 
Solid curve: £ = 0.5 nm; and dotted curve: £ = 0.8 nm. 



Vcnp = 20 V (60 V), respectively, and the solid (dashed) 
curves are from our calculation with Vcnp = 20 V 
(60 V). The fitting parameters are n° = 0.5 x 10 12 cm~ 2 , 
Ag = 0.052 meV and 6 = 2 meV. The variation of £ with 
T is shown in Fig. [TJ^c), where the solid (dashed) curve 
is for Vcnp = 20 V (60 V). The variation of nf with T is 
also shown in Fig. [TJc) with the scale on the right-hand 
side of the frame, where the dotted (chain) curve is for 
Vcnp = 20 V (60 V). In this fitting, the DP mechanism 
is dominant and the spin-flip scattering due to the RRF 
can be neglected. In fact, the calculation with similar pa- 
rameters in Fig. [5] has indicated that the spin relaxation 
time caused by the spin-flip scattering due to the RRF 
is very long. With the same parameters, we further fit 
the dependence of spin relaxation time on diffusion coef- 
ficient at 100 K in Fig. [S] In consistence to the fittings in 
the previous section, the correlation length of the RRF is 
also found to decrease with increasing electron density. It 
is noted that the open squares in the figure are the data 
measured for holes or near the charge neutrality point 26 
and hence are not considered in our fitting. 



D. A nonmonotonic dependence of t s on D 

In the experiments of both Riverside 26 and 
Groningen 2 ^ groups, the spin relaxation time t s is 
observed to depend on the diffusion coefficient D 
monotonically. However, very recently, a nonmonotonic 
dependence of r s on D with the increase of carrier 
density at T = 4.2 K was reported by Jo et al.. 27 



9 



Exp.: V CNP =20 V ■ 
V CNP =60V O 



nf=0.5xl0 12 cm" 2 




50 100 150 200 250 300 

T(K) 



FIG. 7: (Color online) Fit to the experimental data of Han 
and Kawakami2& with the single-sided adatom model, (a) and 
(b): the temperature dependence of t 3 and D. The squares 
(open circles) are the experimental data with Vcnp = 20 V 
(60 V), and the solid (dashed) curve is from fitting with 
Vcnp = 20 V (60 V). (c): the temperature dependences of 
£ and nf (on the right-hand side of the frame). The solid and 
dotted curves are for the case with Vcnp = 20 V, and the 
dashed and chain curves are for the case with Vcnp = 60 V. 



Although this phenomenon is reported at the hole band, 
we can treat it at the electron band by our model due 
to the electron-hole symmetry of band structure in 
graphene. In Fig. [91 we fit the experimental data by the 
single-sided adatom model with nf = 0.5 x I0 12 enr 2 , 
Aq = 0.127 meV and 5 = 4 meV. The closed squares 
(solid triangles) are the experimental (fitting) data of 
the spin relaxation time r s , and the open squares (solid 
circles) are the experimental (fitting) data of the diffu- 
sion coefficient D (with the scale on the right-hand side 
of the frame). The inset shows the density dependence 
of £ (open triangles with the scale on the right-hand 
side of the frame) and nf (solid circles). Due to the 
slower decrease of £ and the faster increase of nf with 
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FIG. 8: (Color online) Fit to the experimental data (the de- 
pendence of spin relaxation time t s on the charge diffusion 
coefficient D with the increase of n e ) of Han and Kawakam: 26 
via the single-sided adatom model. The solid squares stand 
for the experimental data and the open circles are from our 
calculation. The inset shows the dependences of £ (with the 
scale on the right-hand side of the frame) and nf on D when n e 
is increased. The open squares are the experimental data for 
holes or near the charge neutrality point, which are not con- 



sidered in our fitting, nf 
and 8 — 2 meV. 



0.5 x 10" cm"', Aq = 0.052 meV 



increasing n e , it is possible for t s to decrease with n e 
when the latter is large enough. 



E. Possible factors affecting the correlation length 
of the RRF 



In our fittings to the experiments, the variation of the 
correlation length £ of the RRF plays an essential role. £ 
is found to decrease with increasing electron density and 
increase with growing temperature. £ is also found to 
decrease with the increase of adatom density when the 
latter is high enough. It is indeed quite possible that £ 
is affected by these factors. For example, the correla- 
tion length might be shortened by the screening of car- 
riers, which is more effective when the carrier density 
is high. It is also possible that with growing tempera- 
ture the adatoms tend to form clusters and enhance the 
inhomogeneity4^ Besides, the puddle size, which mea- 
sures the correlation length of the short-range random 
potential in graphene, decreases with increasing density 
of the charged impurities^ The similar feature is also 
expected in the current study, i.e., £ decreases with in- 
creasing adatom density when the latter is large enough. 
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FIG. 9: (Color online) Fit to the experimental data (the de- 
pendence of spin relaxation time t s and charge diffusion coef- 
ficient D on electron density n e ) of Jo et al^ via the single- 
sided adatom model. The closed squares (solid triangles) are 
the experimental (fitting) data of the spin relaxation time t s , 
and the open squares (solid circles) are the experimental (fit- 
ting) data of the diffusion coefficient D (with the scale on the 
right-hand side of the frame). The inset shows the density 
dependence of £ (with the scale on the right-hand side of the 
frame) and nf. nf = 0.5 x 10 12 cm -2 , Af, = 0.18 meV and 
6 = 6 meV. 



DISCUSSION AND CONCLUSION 



Discussion on the possible dominant spin 
relaxation mechanism 



We summarize our numerical fittings to the experi- 
ments and discuss the possible dominant spin relaxation 
mechanism. By fitting the DP-like behaviour with the 
increase of adatom density observed by the Riverside 
groups we find that the DP mechanism is the dominant 
one, only with which can the experimental phenomenon 
be understood. It is noted that the correlation length £ 
of the RRF is supposed to be constant at the low density 
regime of adatoms, and when £ is fixed F^p always in- 
creases with the adatom density. Therefore, other kinds 
of attempts with the spin-flip scattering due to the RRF 
being dominant fail to reproduce the experimental phe- 
nomenon and can be ruled out. However, the EY-like 
behaviour with the increase of electron density observed 
by the Groningen group^ 3 . can be fitted by our model with 
either the DP mechanism or the spin-flip scattering due 
to the RRF being dominant when the decrease of £ with 
increasing electron density is considered. Nevertheless, 
the fact that Riverside group has also observed the simi- 
lar EY-like behaviour in their samples very recently^ in 



dependence of the spin relaxation^ suggests that the 
DP mechanism is dominant, but exhibiting EY-like prop- 
erties with the increase of electron density. The similar 
experimental phenomenon on the electron density depen- 
dence of the spin relaxation from the Groningen 2 ^ and 
Riverside 26 groups further suggests that the DP mech- 
anism also dominates the spin relaxation in the experi- 
ment of Groningen groupj^ 3 . Consequently, with the DP 
mechanism being dominant in graphene, the RRF leads 
to spin relaxation which exhibits either DP- or EY-like 
properties in the experiments. 



B. Conclusion 

In conclusion, we have studied electron spin relaxation 
in graphene with random Rashba field by means of the 
kinetic spin Bloch equations, aiming to understand the 
main spin relaxation mechanism in graphene of the cur- 
rent experiments. Different from the previously studied 
case by Zhou and Wu where no adatoms are considered 
and the mobility is relatively highjii^I the electron mobil- 
ity investigated in the present work is at least one order of 
magnitude smaller due to the extrinsic factors caused by 
the ferromagnetic electrodes used in the spin relaxation 
measurements^ -23 ! 26 ' 27 e.g., the adatoms. We set up a 
random Rashba model to incorporate the contribution 
from both the adatoms and substrate. In this model, 
the charged adatoms on one hand enhance the Rashba 
spin-orbit coupling locally and on the other hand serve 
as Coulomb potential scatterers. 

Based on the random Rashba model, the analytical 
study on spin relaxation in graphene is performed. The 
average of the random Rashba field leads to spin relax- 
ation limited by the D'yakonov-Perel' mechanism, which 
is absent in the work by Dugaev et al.^ while the ran- 
domness of the random Rashba field causes spin relax- 
ation by spin-flip scattering, serving as an Elliott- Yafet- 
like mechanism. With the increase of adatom density, the 
spin relaxation caused by the spin-flip scattering due to 
the random Rashba field always shows an Elliott- Yafet - 
like behaviour, whereas the D'yakonov-Perel' mechanism 
can exhibit either Elliott- Yafet- or D'yakonov-Perel'-like 
one. When all the other parameters are fixed, with the 
increase of electron density the spin relaxation rates due 
to both mechanisms increase; Nevertheless, the spin re- 
laxation rate determined by the spin-flip scattering due 
to the random Rashba field is insensitive to the temper- 
ature whereas that determined by the D'yakonov-Perel' 
mechanism becomes insensitive to the temperature when 
the electron-impurity scattering is dominant. However, 
both mechanisms are sensitive to the correlation length 
of the random Rashba field, which may be affected by 
the environmental factors such as electron density and 
temperature. 

We further carry out numerical calculations and fit the 



experiments of Riversid o 24 ' 26 and Groningen 2 ^ groups, 
combination with the observation of the adatom density which show either D'yakonov-Perel'- or Elliott- Yafet 
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like property. By fitting and comparing these experi- 
ments, we suggest that the D'yakonov-Perel' mechanism 
dominates the spin relaxation in graphene. With the 
D'yakonov-Perel' mechanism being dominant, the ran- 
dom Rashba field leads to spin relaxation which exhibits 
either D'yakonov-Perel'- or Elliott- Yafet-like properties. 
Besides, a latest reported nonmonotonic dependence of 
t s on D by Jo et al^ is also fitted by our model with 
the D'yakonov-Perel' mechanism being dominant. 
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Appendix A: Spin-flip scattering terms 

The spin-flip scattering terms axeM- 
dtP^(t)\ f scat = ~Y1 l A k-k'| 2 <5(£Mk-^k') 

X [/^kW^kk'Vk'k + Vkk'Vk'kP M k(*) 

-2Vkk'P M k'(*)Vk'k] 



2tt 



|Ak-k'| 2( 5( e Mk — e ^k') 



x [p M k(t) - Vkk'P M k'(*)Vk'k], (Al) 



where 



[A(r) - A ]e- iq r dr, 



(A2) 



satisfying A q= o = and A q ^o = A q . When the mean free 
path I is much larger than the correlation length £ of the 
fluctuating Rashba field, |A q | 2 can be approximated by 
its statistical average as follows, 



|A q | 2 



x e 



drdr'([A(r)-A ][A(r')-A ]> 

iq- (r— r') 



drdr'C{r - r')e- iq '( r - r '' = C q . (A3) 
Eqs. and JX3| lead to Eq. iJIDj). 



Appendix B: Analytical study of spin relaxation in 
graphene 

We present the analytical study of spin relaxation in 
graphene in detail. By expanding Sk(t) as Sk(£) = 
J2i S k( t ) elWk ' one obtains from Eq. (JT^J) the following 
equations, 



dtS l kx {t) = t E sl kV°(t) - (4- + i)sLW 
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Here |(7 q | 2 and C q depend only on |q| = 2fcsin|. It is 
noted that \ = -4r and -t- = -4r . It is also noted that 

< V T L r ks l 

■T7 = and Tfc is in fact the momentum relaxation time 
T p (k) limited by the electron-impurity scattering. 

By retaining the lowest three orders of S k (t) (i.e., terms 
with I = 0, ±1) in Eqs. (lBl ) -(lB3" l) . one obtains 
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where 
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(B7) 
(B8) 
(B9) 

(BIO) 



As the spin flipping rate 1/t? is much smaller than the 
momentum relaxation rate 1/TjJ (in graphene 1 /rl is usu- 
ally of the order of 10 ps _1 ; even if l/rjj reaches the 
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experimental value ~ 0.01 ps -1 , t^/t° s is still as small 
as 10~ 3 ) and SJr 1 are smaller terms compared to S° in 
the strong scattering limit, we approximate F and Q 
as F w — -p- and Q « 0. With initial conditions, e.g., 

S' fc (0) = i,0,S fc ° z (0)) T (I = 0,±1), Eq. JBi can be 
solved as [those of S^ 1 ^) are away from our interest and 
are not shown here] 
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In the strong scattering limit with c z -C 1, <Sfc z (i) ~ 
SL(0)e- rJt where T+ a 2/r° s + iX^/h 2 . Conse- 
quently, one obtains 



S° k (t) = e-^(0,0,S o kz (0)f 



with 



(B14) 



(B15) 



Similarly, with S^O) = S w (S° kx (0), 0, 0) T or S' fe (0) = 
S lo (0,S° ky (0),0) T , Sl{t) is solved to be 

S°(<) = e - r ^(5L(0),0,0) T , (B16) 



and 



S° fc (i)= e - r ^(0,^(0),0) T , 
respectively, with r x = Y y = T z /2. 



(B17) 
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